Mon. Not. R. Astron. Soc. 000,[T][T2](2012) Printed 10 April 2013 (MN BTfiX style file v2.2) 



Propeller-driven outflows from an MRI disc 



Patrick S. Lii 1 *, Marina M. Romanova 1 , Galina V. Ustyugova 2 , Alexander V. Koldoba 2 , 
Richard V.E. Lovelace 1 ' 3 



o 

(N 

< 

0\ 



6 

in 



> 

o 

(N 

O 
m 

> 

• i— I 

X 

$H 



1 Department of Astronomy, Cornell University, Ithaca, NY 14853-6801 

2 Keldysh Institute for Applied Mathematics, Moscow, Russia 

3 Department of Applied and Engineering Physics, Cornell University, Ithaca, NY 14853-680 



10 April 2013 



ABSTRACT 

We present the results of axisymmetric simulations of MRI-driven accretion onto a rapidly ro- 
tating, magnetized star accreting in the propeller regime. The stellar magnetosphere corotates 
with the star, forming a centrifugal barrier at the disc-magnetosphere boundary which inhibits 
matter accretion onto the star. Instead, the disc matter accumulates at the disc-magnetosphere 
interface and slowly diffuses into the inner magnetosphere where it picks up angular mo- 
mentum and is quickly ejected from the system as an outflow. Due to the interaction of the 
matter with the magnetosphere, this wind is discontinuous and is launched as discrete plas- 
moids. If the ejection rate is lower than the disc accretion rate, the matter accumulates at the 
disc-magnetosphere boundary faster than it can be ejected. In this case, accretion onto the star 
proceeds through the episodic accretion instability in which episodes of matter accumulation 
are followed by simultaneous accretion and ejection. During the accretion phase of this insta- 
bility in which matter flows onto the star in funnel streams, we observe a corresponding rise 
in the outflow rate. Both the accretion and ejection processes observed in our simulations are 
highly non-stationary. The stars undergo strong spin-down due to the coupling of the stellar 
field with the disc and corona and we measure the spin-down timescales of around 1 Myr for 
a typical CTTS in the propeller regime. 
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1 INTRODUCTION 

Young, newly formed protostars rotate rapidly and may be in the 
"propeller" regime of disc accretion in which the angular veloc- 
ity of the stellar magnetosphere is larger than the angular veloc- 
ity of the inner disc (Illarionov & Sunyaev 1975). Examples of 
other propelling stars include accreting fast-rotating neutron stars 
as well as white dwarfs in cataclysmic variables (Stella et al .|1986| 
|Cui|1997||Alpar|2001[|Patruno et al.|2009| l. In these stars, accre- 
tion in the inner disc is determined by the relative locations of 
the disc-magnetosphere boundar}r|r m ag and the corotation radiusrl 
?"cor — {GM^/nt) 1 ' 3 . If r mag <r cor , then the magnetosphere ro- 
tates slower than the inner disc and the disc is disrupted, causing 
the matter to accrete onto the star as a funnel flow (Lamb et al. 
|1973[|Elsner & Lamb|1977||Romanova et al.|2002) . In the oppo- 
site situation where r mag >r cor , the magnetosphere rotates more 
rapidly than the inner disc. In this case, the matter in the inner disc 



* E-mail of corresponding author: pslii@astro . Cornell . edu 

1 Determined by the balance between the magnetic and matter pres- 
sures in the inner disc, so that the kinetic plasma parameter at the disc- 
magnetosphere boundary fi\ = (p + pv 2 )/(B 2 /8n) ~ 1. 

2 The radius at which the angular velocity of the star f!„ equals the local 
Keplerian angular velocity. 



may acquire angular momentum from the magnetosphere, leading 
to super-Keplerian rotation of the matter which may in turn drive 
an outflow ( jLovelace et al.| 1999). This is known as the propeller 
regime of accretion. 

The propeller regime has been studied analytically using mag- 
netohydrodynamics (MHD, e.g. Davies et al.|1979| |Li & Wickra 



masin ghe|T99"7l |Lovelace et al.|l999||I1disanov||2002||Rapp'aport' 

et al. 2004 , Eksi et al. 2005 ) as well as numerically in MHD sim 



ulations (Wang & Robertson 1985, Romanova et al. 2004, 2005; 
Ustyugova et al. 2006, Romanova et al. 2009). Numerical stud- 



ies by [Romanov a et aTj l 2005) and Ustyugova et al.| pOOo} con- 
firmed that propelling stars can launch strong outflows as predicted 
by [Lovelace et af] ( |1999| >. The outflows are found to have two 
main components: (1) a conical-shaped, centrifugally-driven mag- 
netospheric wind into which much of the disc matter flows and 
(2) a high-velocity, magnetically driven and collimated axial jet 
into which a significant amount of energy and angular momentum 
flows. These previous simulations incorporated an a-prescription 
(Shakura & Sunyaev 1973) to model the viscosity and magnetic dif- 
fusivity of the disc plasma. Additionally, the simulations were per- 
formed on a spherical grid assuming top-bottom symmetry across 
the equatorial plane as well as axisymmetry about the rotation axis. 
More recently, one-sided outflows and opposite sided funnel flows 
were found in axisymmetric simulations of accretion of a viscous 
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and diffusive a-disc onto a rapidly rotating magnetized star for 
conditions where symmetry about the equatorial plane was not im- 
posed ( |Lovelace et al.| 2010). 

Here, we study the problem of propeller outflows utilizing a 
new code which has a high spatial resolution in the region occupied 
by the disc. The higher resolution grid resolves the MHD turbu- 
lence generated by the magnetorotational instability (MRI, Balbus 
|& Hawley|1 991). The MRI-driven turbulence causes accretion of 
the disc and for this reason we do not include either a viscosity or 
diffusivity coefficients, considering only an ideal MHD disc. Ear- 
lier, this code was used to study MRI-driven accretion onto slowly 
rotating stars (Romanova e t al.|20lT) . 

MRI-driven accretion has been studied in shearing boxes as 
well as in global simulations onto non-magnetized stars or black 
holes (e.g.,|Hawley et al.[T995l|Stone et al.|1996 . Armitag ep998l 
Gammie & Men ou|1998[|Haw ley 2000[ |Hawley et al.|2001||stone| 
& Pringle|200Tp iawley & Kroli kl2002[|Beckwith et al.|2009> . Re- 



cently, our group was able to investigate MRI-driven accretion onto 



a magnetized star, first under axisymmetric conditions (Romanova 



et al.pOll) and later in global 3D simulations (Romanova et al. 



2012 1. It was shown that the disc matter flows onto the star through 



funnel streams (as in |Romanova et aLp 002 2003, |2004| l. 

In the present work, we model and study MRI-driven disc ac- 
cretion onto a star in the propeller regime where r mag > r cor . In 
this regime we observe episodic funnel accretion which is accom- 
panied by centrifugally-driven outflows launched from the disc- 
magnetosphere boundary. The aim of this work is to analyze the 
accretion and outflow processes in the propeller regime. In O we 
discuss the numerical model and initial conditions of the simula- 
tions. In ipl we present an overview of the three models which we 
consider in this work; in §3.1| we investigate the launching mech- 
anisms; in j ]3.3| we measure the spin-down rates for each of our 
models; and in ij4] we discuss the structure of the outflows and the 
application of our models to real systems. 



2 THE NUMERICAL MODEL 

Code description: We use a Godunov-type numerical method 
with a five-wave Harten, Lax, and van Leer (HLL) Riemann solver 
similar to the HLLD solver developed by Miyoshi & Kusano 
(2005). The MHD variables are calculated in four states bounded 
by five MHD discontinuities: the contact discontinuity, two Alfven 
waves and two fast magnetosonic waves. Unlike Miyoshi & Ku- 
|sano| (2005), our method solves an equation for the entropy in- 
stead of the full energy equation. This approximation is valid in 
cases such as ours where strong shocks are not present. In order to 
avoid non-physical solutions of the Riemann problem, we perform 
a correction procedure on the fast magnetosonic wave velocities 
in order to maintain the gap between these waves and the Alfven 
waves which propagate behind the fast magnetosonic waves. We 
ensure that the magnetic fields are divergence-free by introducing 
the (^-component of the magnetic field potential which is calcu- 
lated using the constrained transport scheme proposed by Gardiner 
& Stone (2005). To minimize the error in the Lorentz force, we 
split the magnetic field into the stellar dipole and calculated com- 
ponents, B — Bdip + B' , omitting terms of the order Bdip which 
do not contribute to the Maxwellian stress tensor (Tanaka 1994). 
No viscosity or diffusivity terms have been included in the MHD 
equations and hence we investigate only accretion driven by the re- 
solved MRI-turbulence. Our code has been extensively tested and 
has been previously utilized to study MRI-turbulent accretion onto 




30 10 20 




20 30 10 20 30 



Figure 1. The Grid. From left to right: (1) grid used in simulations, (2) 
same grid but with l/10th of the grid lines shown for clarity; (3) initial 
distribution of density (background) and magnetic flux (lines); (4) same but 
for time t = 500. 



a magnetized star (see tests and other details in Romanov a et al.| 
|20TT] >. 

Grid description: The axisymmetric grid uses cylindrical (r, z) 
coordinates with mesh compression towards the disc and towards 
the «-axis such that there are a larger number of cells in the disc 
plane and near the star (see Fig. IT] left panel). In the models pre- 
sented here, we use a non-uniform grid with dimension 250 x 432 
cells corresponding to a grid that is 33 by 70 stellar radii in size 
(0.31 AU by 0.65 AU if P» = 2B Q ). At r = 20, the number of 
grid cells which cover the disc in the vertical direction is about 200. 

Reference units: The simulations are performed in dimensionless 
units and are applicable to stars over a wide range of scales so long 
as the magnetospheric radius r mag is not very large compared to 
the radius of the star P* (i.e., r mag = (4 — 6)7?,). We choose the 
values of the stellar mass Af„, radius P,, magnetic field B* and 
dimensionless magnetic moment p, = 10 and derive reference val- 
ues using the initial values as a base. The derived reference units 
are as follows: length Ro = B„, magnetic moment fio — BoBo, 
magnetic field Bo = B*/p, with the dimensionless magnetic mo- 
ment jl = 10 for our simulations, velocity no = \/GM,/Rq 
(the Keplerian orbital velocity at the surface of the star), time pe- 
riod Po = 2ttRo/vo (the Keplerian orbital period at the surface 
of the star), angular velocity fio = vo/Ro, force per unit mass 
/o = Vq/Ro, pressure po = B , density po = Po/vq, temperature 
To = Po/ po, accretion rate Mo = poVoRo an d angular momen- 
tum flux Lo = MovoRo- Table 111 shows sample reference values 
for three different types of accreting stars: to apply the simulation 
results to a particular class of star, multiply the dimensionless value 
by the reference value. 

For the remainder of this work, all values and variables are 
given in terms of the normalized units, except where they have ex- 
plicitly assigned physical units. In all following plots, time and dis- 
tance are always given in units of Po and Po. 

Initial conditions: The initial conditions for the hydrodynamic 
variables are similar to those taken in our previous works (e.g. Ro- 
|manova et a l. 2002 2011) where the initial density and entropy 
distributions have been calculated from the force-balance of differ- 
ent forces which include gravity, centrifugal and pressure forces. 
The disc is dense and cold while the corona is about 1000 times 
hotter and of lower density. 

We study three main cases of propelling stars in different rota- 
tional regimes with corotation radii r COI = 1.3, 1.5, and 2 (see Table 
[2] for other parameters). At the start, the inner edge of the disc is 
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CTTS 




White Dwarf 


Neutron Star 


initial 


M, [M ] 


0.8 




1 


1.4 


R* 


2R© 




5000 km 


10 km 


B.[G] 


1000 




5 x 10 5 


5 x IO 8 


derived 


Ho [cm] 


1.39 x 


10 11 


5 x 10 8 


1 x IO 6 


Bo[G] 


100 




5 x IO 4 


5 x IO 7 


no [cm s ] 


2.76 x 


10' 


5.15 x 10 8 


1.36 x IO 10 


Po 


0.37 d 




6.1s 


0.46 ms 


n [s- 1 ] 


1.99 x 


10~ 4 


1.03 


1.36 x IO 4 


/oCdyg" 1 ] 


5.49 x 


10 3 


5.31 x 10 8 


1.86 x IO 14 


po [dy cm -2 ] 


1.0 x 10 4 


2.5 x 10 9 


2.5 x IO 15 


po [g cm -3 ] 


1.31 x 


10" 11 


9.42 x 10~ 9 


1.35 x 10~ 5 


/in [ergG -1 ] 


2.69 x 


10 35 


6.25 x IO 30 


5 x IO 25 


M [M© yr" 1 ] 


1.11 X 


io- 7 


1.92 x IO" 8 


2.91 x IO" 9 


i [g cm 2 s" 2 ] 


2.69 x 


10 37 


3.13 x IO 35 


2.5 x IO 33 



Table 1. Reference units. Calculated reference values for three different 
types of accreting stars. We choose the values of the stellar mass M * , radius 
R t , magnetic field B*, and dimensionless magnetic moment /J = 10 and 
derive the other reference values from these initial values. 



Top and bottom boundaries: all variables have free boundary con- 
ditions along the top and bottom boundaries. In addition, we imple- 
ment outflow boundary conditions to prohibit matter from flowing 
back into the simulation region once it leaves. 
Outer side boundary: the side boundary is divided into a "disc re- 
gion" (|zj < Zdisc) and a "coronal region" (\z\ > ZdiscX w ' tn 



Zdi, 



h(R ut) 



GM„ 



<M#o 



E 



Rl 



where R ou t is the external simulation radius. The matter along the 
disc boundary (|z| < Zdisc) I s allowed to flow in with a small radial 
velocity 



— j| 



P 



2 pvk(Rou 



5 = 0.02, 



and with a poloidal magnetic field corresponding to the calculated 
magnetic field at r = R out . The remainder of the variables are 
allowed to have free boundary conditions. The coronal boundary 
(\z\ > Zdisc) has the same boundary conditions as the top and bot- 
tom boundaries. 



'"cor 


Q* 


P* 


A 


1.3 


0.67 


9.31 


10 


1.5 


0.54 


11.5 


10 


2.0 


0.35 


17.8 


10 



Table 2. Models and model parameters. In dimensionless units, r cor is 
the corotation radius, £2, is the angular velocity of the star, P» is the period 
of stellar rotation and fi is the dimensionless magnetic moment. 



placed at r — 10 and the star rotates slowly with £2^=0.032 (corre- 
sponding to r COI = 10) such that the magnetosphere and inner disc 
are initially corotating with the same angular velocity. This condi- 
tion helps to ensure that the magnetosphere and disc are initially in 
a near-equilibrium at the disc-magnetosphere boundary. From to 
100Po, the star is gradually spun up from Qj to the final rapidly- 
rotating state with angular velocity fi* (given in Table[2]l. 

Initial magnetic field configuration: Initially, the disc is threaded 
with the dipole magnetic field of the star. We also add a small "ta- 
pered" poloidal field inside the disc (see third panel in Fig.fT) which 
is given by 



,T, BoT f Z 

W = COS TV 

2 V 2h 



h = 



GAL 
$ c (r) - E 



where h is half-thickness of the disc, & c {r) = kGAh/r where k 
is a Keplerian parameters and £ is a constant of integration in the 
initial equilibrium equation (see Romano va et al. 2002, 201 ll. The 
magnetic field in the disc has the same polarity as the stellar field 
at the disc-magnetosphere boundary. 

Boundary Conditions: Stellar surface: all variables on the sur- 
face of the star have "free" boundary conditions such that 
d(...)/dn — along the entire surface. We prohibit the outflow 
of matter from the star (i.e. we do not allow for stellar winds) and 
adjust the matter velocity vectors to be parallel to the magnetic field 
vectors. This models the frozen-in condition on the star. 



3 We take k slightly greater than one to balance the disc pressure gradient 
(k=l+0.02). 



3 RESULTS : OUTFLOWS IN THE PROPELLER REGIME 

We study three main cases of stars accreting in the propeller regime 
with different rotational periods corresponding to r COI = 1.3, 1.5 
and 2. For a representative T Tauri star (M» = 0.8 Mq, R* = 
2 Rq) these corotation radii imply rotational periods of 3.4, 4.3, 
and 6.6 days, respectively. The dimensionless parameters applica- 
ble to other types of stars can be found in Table [2] For stars ac- 
creting in the propeller regime, several processes are expected to 
occur: firstly, matter which diffuses into the rapidly rotating mag- 
netosphere from the disc will be rapidly accelerated past the local 
Keplerian velocity and driven out as a propeller wind. In the re- 
gions above the star where the matter density is low, the magnetic 
field lines can easily inflate into the corona and the star may in- 
stead drive a magnetically dominated Poynting jet. Secondly, the 
rapidly rotating star should quickly spin down as angular momen- 
tum flows out from the star through the magnetic field lines con- 
necting it to the corona and disc. In this section, we study both 
of these effects: in j ]3.1| we investigate the matter-dominated wind 
driven from the disc by the propelling star, in |3.2| we discuss the 
magnetically-dominated Poynting component of the outflow, and 
in j ]3.3| we measure the stellar spin down rate associated with the 
outflow in each of the models. 



3.1 Outflows in the propeller regime 

The propeller outflow consists of two main components: a matter- 
dominated propeller wind launched from the disc-magnetosphere 
boundary and a magnetically-dominated Poynting jet which flows 
along the field lines connecting the star to the slower rotating 
corona. In this section, we discuss the matter dominated wind com- 
ponent. 

In order for the propelling star to eject matter from the disc, 
the magnetospheric velocity must exceed the Keplerian velocity 
at the disc-magnetosphere boundary, f(r mag ) > Wis-(r mag ). If 
this condition is satisfied, the rapidly rotating magnetosphere acts 
as a centrifugal barrier and the matter in the disc is prevented 
from accreting onto the star. Instead, the matter accumulates at the 
disc-magnetosphere boundary and diffuses into the magnetosphere 
where it picks up angular momentum and is ejected as an outflow. 
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Figure 2. Mass fluxes in the three models. The red lines show the accretion rate onto the star, M* , and the dark blue lines show the outflow rates M w i n d . 
For reference, the time averaged outflow rates (M w ; n< j) are also plotted as black dashes. 



The matter diffusion therefore determines the outflow rate and in 
Appendix [Al we further investigate the role of the diffusivity on 
the outflow rate in a special set of non-ideal MHD simulations. In 
the models, the outflow launched in this non-accreting state is non- 
stationary but proceeds continuously: when the matter is ejected, 
it inflates the field lines in the outer magnetosphere which expand 
and then quickly reconnect. This cyclic expansion-reconnection re- 
curs on the dynamical timescale of the inner disc and the matter is 
ejected as discrete plasmoids similar to the clumpy winds observed 
in previous studies of propeller outflows (e.g Romanova et al. 2005 
|Ustyugova et al.|2006| >. 

If the diffusion rate is low and the disc's accretion rate ex- 
ceeds the outflow rate, the accreting matter is halted at the disc- 
magnetosphere boundary. Eventually, this accumulated matter ac- 
cretes onto the propelling star through the episodic disc instabil- 
ity first described by Spruit & Taam ( 1993) and studied further by 
D' Angelo & Spruit (2010). The basic mechanism behind the insta- 
bility is as follows: 

(i) The matter from the disc is blocked from accreting by 
the centrifugal barrier and accumulates at the disc-magnetosphere 
boundary, gradually compressing the stellar magnetosphere to- 
wards the star. 

(ii) The matter continues to compress the magnetosphere toward 
the star until the gravitational acceleration is larger than the cen- 
trifugal acceleration causing the matter to accrete onto the star as a 
funnel flow. 

(iii) With the reservoir of matter depleted, the magnetosphere 
quickly re-expands out and shuts off the accretion, restarting the 
cycle anew. 

Episodic accretion results in a "spiky" accretion rate in which any 
accretion onto the star occurs in quasi-cyclic bursts. The accretion 
rate M» in the top panels of Fig.blshow these characteristic bursts 
of matter in each of the three propeller models which we consider 
here. We measure the wind outflow rate, A/ W i n d, by calculating the 
matter flux through a cylindrical surface centered on the star with 
radius r = 20 and height z = 40 and only measuring the fast 
component of the flow with v p > 0.2. Figure [2] shows that the 
outflows (Af w i n( j) are highly non-stationary and vary strongly from 
model to model. In the simulations, these bursts of accretion tend to 
be accompanied by a burst of outflowing matter as well, as shown 
byM wind . 

To get a clearer picture of the accretion and outflow rates, we 
calculate accretion rate time averages 



(M{t)) = 






[2] the fastest rotating star has the largest overall outflow rate while 
the slowest rotating star has the smallest. We calculate accretion to 
ejection ratios <M*)/(M W md) of 0.144, 0.345, and 0.964 for the 
r cor = 1.3, 1.5, 2 cases respectively. These values indicate that the 
faster rotating stars in the r cor = 1.3 and 1.5 cases eject most of the 
accreting matter into an outflow before it ever reaches the surface of 
the star. In contrast, the weaker r cor = 2 propeller does not launch 
the matter as efficiently and (M*) ~ (M w i n d) for this case. 



(1) 



The time averaged (M w i n d) is shown as a dashed line in Figure 



Quasi-steady outflows in the non-accreting phase 

During the accumulation phase of the accretion-ejection cycle, 
some of the matter at the disk-magnetosphere boundary diffuses 
through the disc-magnetosphere boundary into the inner magneto- 
sphere where it picks up angular momentum and is ejected as a 
centrifugally-driven magnetospheric outflow. Figure [3] shows three 
time snapshots of these centrifugally launched outflows in the r cor = 
1 .3 model. The top panels show the logarithmic poloidal matter flux 
p|v p |, poloidal velocity vectors v p , and the magnetic field lines 
near the star. In order for the matter to be driven into an outflow, it 
must overcome the local gravitational force as well as the magnetic 
tension in the magnetospheric field lines. The matter which flows 
out from the magnetosphere opens up the field lines, inflating the 
magnetosphere in the direction of the outflow (t = 310). As the ac- 
cumulated matter flows away, the magnetic tension causes the field 
lines to reconnect (t = 312) and truncate the outflow into discrete 
lumps of matter which flow away at roughly 1.5 times the local 
escape velocity. As in previous studies of propeller outflows jRo-| 
manova et al.|2005^ , we observe a half-opening angle of roughly 
45° at the base of the propeller wind with evidence of gradual col- 
limation at larger distances. The bottom panels of Figure [3] show 
snapshots of the angular velocity £7 alongside the poloidal velocity 
vectors v p and stellar magnetic field lines. The large contrast in the 
angular velocities of the magnetosphere and inner disc is evident in 
the figure. Additionally, the reservoir of matter which builds up at 
the disc-magnetosphere boundary can be seen in the geometry of 
the magnetic field lines. 

An interesting feature of these propeller-driven outflows is that 
they are asymmetric. Previous studies of outflows from propelling 
stars have always assumed equatorial symmetry. Here, we see that 
one-sided centrifugally-driven outflows are launched even if the 
magnetosphere is initially top-bottom symmetric, similar to the 
one-sided magnetic outflows observed in simulations by |Lovelace| 
[eTaLlpOTO) . 
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Figure 3. The non-accreting outflow mode. Top panels: three snapshots of the r cor =1.3 simulation region showing the plasmoids launched during the non- 
accreting phase of the episodic accretion cycle. Logarithmic maps of pv p axe plotted with field lines and velocity vectors superimposed. Middle panels: the 
angular frequency Q of the matter at those same three moments in time. Bottom panels: matter and angular momentum fluxes onto the star and into the wind. 
The vertical dotted lines denote times corresponding to the three snapshots shown. 



3.1.1 Outflows driven in the accreting phase 

In addition to the outflows in the non-accreting stage, the simula- 
tions show a stronger transient outflow associated with episodes of 
funnel flow accretion from the disc onto the star. Figure H] shows 
four snapshots of one such accretion event in the r cor = 1.3 model. 
At t = 35 1, the matter reservoir at the disc-magnetosphere interface 
has compressed the magnetosphere inward enough for accretion to 
proceed onto the star. By t = 356, a funnel flow has formed and the 
accretion and outflow rates have peaked: as the disc matter flows 
along the field lines toward the star, the contact area between the 
matter and magnetosphere is extended, permitting matter to diffuse 
through into the inner magnetosphere. There, the matter quickly 
picks up angular momentum from the field and is driven out into an 
outflow. Again, the diffusion of the matter across field lines plays 
a crucial role in loading the outflow (see Appendix [A] for further 
discussion). Once the reservoir of matter compressing the magne- 
tosphere is depleted by accretion onto the star, the magnetosphere 
re-expands, halting the funnel flow of matter and reverting the out- 
flow back to the initial weaker state. This episodic accretion event 
creates a short burst in the outflow rate occurring on a few dynam- 



ical timescales of the very inner disc. The stronger outflows in this 
phase tend to be more collimated than in the weaker, non-accreting 
phase with narrower half-opening angles ranging from 20-40° at 
the simulation boundary, depending on the amount of matter ac- 
creted onto the star. 

In the simulations presented here, the disc accretion rate is 
high enough to compress the magnetosphere and drive the episodic 
accretion instability which leads to the spikes in the outflow rate. 
In stars with lower accretion rates or stronger fields, the magneto- 
sphere may be very large compared to the disc thickness and ac- 
cretion onto the star may be completely blocked. In these cases, 
the stronger outflows associated with the accretion may not be ob- 
served at all as funnel flows onto the star are prohibited. 



3.2 Poynting flux outflow 

In addition to the matter-dominated propeller wind from the disc- 
magnetosphere boundary, we observe a magnetically dominated 
Poynting jet which flows along the helically wound field lines con- 
necting the star with the slower rotating corona. In the jet, mag- 
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Figure 4. The accreting outflow mode. Top panels: four time snapshots of the propeller outflow during an accretion event. A large amount of matter is ejected 
into the outflow when accumulated matter funnel flows around the magnetosphere from the disc onto the star. Middle panels: the angular frequency Q, velocity 
vectors and magnetic field lines at the same four moments in time. Bottom panels: the mass and angular momentum fluxes onto the star and into the wind. The 
vertical dashed lines correspond to the four time snapshots. 




Figure 5. Energy flux densities. Left panel: energy flux density carried by 
the magnetic fields. Middle panel: energy flux density carried by the mat- 
ter. Right panel: logarithmic matter fluxes, poloidal field lines, and velocity 
vectors. All plots show the r cor =1.3 model at t=400 during the accreting 
phase of the cycle. 



netic pressure accelerates the matter away from the star while the 
magnetic hoop stress simultaneously collimates it (Ustyugova et al. 
|2000| [Lovelace et al. 2002). The jet is strongest during the "accret- 
ing" phase when the wind outflow is also the strongest because both 
phenomena are associated with the inflation and ejection of mag- 



netic flux from the stellar magnetosphere. The magnetic Poynting 
outflow is well collimated with half-opening angles of 5-15° at the 
edge of the simulation region and the inner, strongly magnetically 
dominated regions of the jet exhibit more collimation than the outer 
regions of the outflow where there is more matter. 

To study the energy outflow in the jet and wind, we measure 
the energy flux densities associated with the field 



'Et 



4tt 
and the matter 

2 

/ nr 



Em — 



E x B 



IBMv, 



(B v)B c 



4tt 



(2) 



pv 



(3) 



7-1 

where E = — v x B/c and 7 = 5/3 for the ideal plasma in 
our models. Figure B] shows the energy fluxes at t = 400 dur- 
ing an episode of enhanced outflow into the wind: the leftmost 
panel shows the energy flux in the magnetic fields, highlighting the 
Poynting outflow from the star. As matter accretes, the enhanced 
magnetic outflow accelerates some of the accreting matter into a 
high velocity (|v p | ~ 2), well-collimated jet. The middle panel 
shows the matter dominated wind also carries away substantial en- 
ergy from the inner disc. The velocity vectors in the rightmost panel 
show that the accreting matter is accelerated to high velocities in 
both the magnetic and matter dominated regions of the flow. Even 
in the matter dominated components, there is a substantial amount 
of magnetic energy flux which helps to partially collimate the wind 
up to 20^10° at the edge of the simulation region. This magnetic 
flux will continue to collimate the flow at even larger distances. 
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Figure 6. Angular momentum flux densities. Top left panel: Angular mo- 
mentum flux density carried by the magnetic fields. Top right panel: Angu- 
lar momentum flux density carried by matter. All plots show the r cor =1.3 
model at t=3 15. 



3.3 Propeller spin down 

The propelling star loses angular momentum through the stellar 
magnetic fields connecting the star to the slower rotating matter in 
the corona and disc. We are able to measure the angular momentum 
lost by the star by integrating the angular momentum flux densities 
over the stellar surface 



L = L m + L f = / dS ■ (F Lm + F 



(4) 



Here the normal vector for the surface dS points inward towards 
the star and FLm and Fu are the angular momentum flux densities 
carried by the matter and magnetic field, given by 



rpv^Vp , 



4tt 



(5) 



Figure [6] shows the angular momentum flux densities in the r cor = 
1.3 model at t = 315. The left panel shows the angular momentum 
being carried by the magnetic fields: a significant part of the an- 
gular momentum flux lies in the Poynting jet indicating that both 
the Poynting jet and propeller wind are important in spinning down 
the star. This spin-down torque is largely mediated by the magnetic 
field lines connecting it to the slower rotating corona and disc. The 
right panel, which shows the angular momentum transport through 
matter, highlights the matter stress induced by MRI in the disc as 
well as the substantial amount of angular momentum carried by the 
matter itself in the propeller-driven wind. 

To measure the long-term spin-down rate of the star, we time 
average the angular momentum fluxes 



(L®) 



s: Lm 



and overplot the averages over the actual fluxes in Figure [7] In 
agreement with the above analysis, these time averages indicate 
that the majority of the stellar angular momentum leaves in the 
outflow or as Poynting flux while the small fraction remaining is 
transported back outward through the disc by MRI. As shown by 



the unaveraged L t , the episodic nature of the accretion results in a 
torque which is episodic in nature as well. 

Using the time averaged angular momentum fluxes, we es- 
timate the spin-down timescale for the stars in each of the three 
propelling stars. In dimensional units, the characteristic spin-down 
time is 



L 



(6) 



where / is the star's moment of inertia. Substituting in L — LLo 
and fi» = fi*f2o (with tildes denoting the dimensionless values) 
and expanding to reference units, we have 
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where we have plugged in values for a 0.8 Mq, 2 Rq CTTS with 
I « 1.23 X 10 55 g cm 2 and a 1 kG surface field (see Table[j]for the 
other reference units). Using this, we can estimate the spin-down 
timescale for each of our cases: 

(i) r cor = 1.3 ((!/») « 2.0), £ sp ,i.3 » 9.65 x 10 5 yr. 
(ii) r cor = 1.5 ((Lfs) « 1.7), V1.5 « 9.10 x 10 5 yr. 
(iii) r COT = 2.0 ((L fa ) « 1.1), i sp , 2 » 9.50 x 10 5 yr. 

The characteristic spin-down timescales are roughly the same for 
all three stars, hovering around one megayear. These characteris- 
tic timescales are roughly constant because the faster rotators have 
larger spin-down rates but also a larger amount of angular momen- 
tum. Table[3]summarizes some of the dynamical properties of each 
of the three models. The derived spin-down timescales are con- 
sistent with measurements of CTTSs rotation periods which show 
that CTTSs older than 1 Myr are already slow rotators. As shown 
in Equation [7] the spin-down times depend on a large number of 
parameters — in particular, the stellar radius _R», stellar magnetic 
field £>», and the dimensionless magnetic moment jl which deter- 
mines the size of the magnetosphere relative to the star. If the young 
star has a larger magnetosphere or radius than the given reference 
star, it will spin down even more rapidly that the 1 Myr timescale 
estimated here. 



4 DISCUSSION 

4.1 Propelling CTTSs 

Newly formed Class T Tauri stars are expected to be in a state 
of rapid rotation and may have strong magnetic fields which in- 
teract with their accretion discs. For this reason, propeller-driven 
outflows are expected in this early stage of stellar evolution. How- 
ever, in this phase, the nascent stars are enshrouded by an optically- 
thick layer of gas and dust, preventing direct observations of the 
inner launching region. Instead, observations of these young stars 
show Herbig-Haro jets emerging from a dense optically-thick pro- 
plyd which hides the star. One of the possible origins of these ob- 
served outflows is the propeller mechanism discussed in this work 
which channels substantial amounts of matter and energy into a 
well-collimated outflow. 
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Figure 7. Angular momentum fluxes in the three models. The dark red and blue lines show the time averaged angular momentum flux onto the star and into 
the wind, respectively. The pink and light blue lines show the corresponding unaveraged fluxes onto the star and into the wind. 
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Table 3. Summary of values for CTTSs in the propeller regime. These values are for a 0.8 Mq, 2i?0CTTS with a 1 kG surface dipole field (see TablefT] 
for other parameters). r C oris the corotation radius, P* is the rotation period of the star, (Af w ; n( j}/((Af w i n( j} + (A/,)) is the outflow ejection efficiency, t Bp 
is the spin-down timescale of the star and S7„ is the stellar spin-down rate. 



In later stages of protostellar evolution (after ~ 10 yrs), the 
thick cloud of gas and dust enshrouding the star clears and the 
young star becomes accessible to direct observation. At this stage, 
measurements show the stars are relatively slowly rotating, with 
periods -~2-10 days despite the fact that the stars continue to grav- 
itationally contract. This is the Classical T Tauri star stage of the 
evolution (CTTSs, or Class II stars). The CTTSs show outflows 
which are much weaker than at the earlier stages of evolution. Pro- 
peller outflows will be present in this stage during episodes of low 
accretion when the stellar magnetospheric radius exceeds the coro- 
tation radius of the star. 

One such example of a propelling CTTS is the well studied 
Class II T Tauri star AA Tau, with inferred mass, radius and age of 
0.7M Q , 2Rq and 1.5 Myr. Spectropolarimetric measurements by 
Donati et al. (2010) indicate that its magnetic field is dominated by 
a 2-3 kG dipole component inclined at 20° to the rotation axis of 
the star. This inclined dipole creates a warp in the inner disc which 
corotates with the star, periodically eclipsing it every 8.22 days, 
indicating that the star has a corotation radius r cor = 7.6-R, (Bou- 
|vier et al.|20 07 1. Based on these parameters and measurements of 
the accretion rate (10 _9 ' 2 AfQ/yr on average over their observation 
epochs), Donati et al. (2010) showed that AA Tau is well within 
the propeller regime of accretion with r mag ~ 2 r cor . As the mag- 
netospheric radius lies well beyond corotation, AA Tau is a prime 
candidate for launching propeller-driven outflows. Imaging studies 
of the AA Tau jet by |Cox et al.| ( |2013) reveal a clumpy high-velocity 
outflow with an unusually wide half-opening angle of 17° within 2" 
of the star. In the same study, a possible counterjet was suggested, 
but not confirmed. 

The microjets observed in young T Tauri stars are known to 
be well collimated at distances greater than 10 AU from the star. 
Observations suggest that they have a layered "onion-skin" struc- 
ture, with a well-collimated, high-velocity central component sur- 
rounded by several less collimated, lower velocity layers. The well- 
collimated central components typically have half-opening angles 
of 2-12° — narrower than the AA Tau jet and narrower than the 
20-40° wind that we observe in our models (e.g., Dougados et al. 
|2000| [Agra-Amboage et al.|20lT| >. In order for the propeller wind 
to explain the jet phenomenon in T Tauri stars, a larger-scale col- 
limation mechanism must be present. The most likely scenario is 



that the propeller wind continues to be self-collimated beyond the 
edge of our simulation region. Despite being matter dominated, 
the propeller wind carries a substantial amount of magnetic en- 
ergy which is responsible for the partial collimation observed in 
our models (see j |3.2) . This magnetic flux will continue to collimate 
the flow at larger distances, resulting in a well-collimated propeller 
outflow far from the star. It may also be possible for the propeller 
wind to be confined by the ram pressure of a slower moving larger- 
scale disc-wind or by the magnetic pressure of large scale poloidal 
fields threading the accretion disc, but it is not clear if these mech- 
anisms are important as the magnetic collimation is already suf- 
ficient prank & Mellema|[T996l |Matt et al.||2003) . A collimated 
propeller outflow may be responsible for the central medium/high- 
velocity components observed in T Tauri jets while the lower ve- 
locity outer layers may be due to a extended magneto-centrifugal 
disc-wind launched at larger distances. Observationally, it is not at 
all clear if the microjet in AA Tau is in fact propeller-driven, but 
the outflow velocity and unusually wide opening angle may be ex- 
plained by a propeller wind which has undergone weak collimation. 



4.2 The one-sided propeller wind 

The geometry of the stellar magnetosphere plays a key role in the 
launching of the propeller outflow. In our axisymmetric models, 
the aligned dipole is initially equatorially symmetric. However, as 
soon as the matter from the disc reaches the magnetosphere, the 
stellar dipole is deflected to one side and this initial top-bottom 
symmetry is broken. This asymmetry helps to drive the one-sided 
outflows observed in this work as it provides a preferential direction 
for the matter accretion and ejection. For our simple dipole magne- 
tosphere, we observe occasional "flip-flops" where the direction of 
the outflow switches sides. Lovelace et al. (2010) studied one-sided 
outflows driven by stars with both dipole and quadrupole compo- 
nents finding that an aligned, purely quadrupolar field could drive 
top-bottom symmetric outflows. However, the addition of even a 
small dipole component to the field caused one-sided outflows to 
be driven from the disc-magnetosphere boundary. As in the simu- 
lations presented here, |Lovelace et aT| |2010) also found that a pure 
dipolar field drives one-sided outflows which can occasionally flip 
sides. 
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The presence of a complex magnetic topology on the star may 
also help to drive one-sided outflows. Most CTTSs have fields 
which are inclined relative to the rotation axis — a configuration 
which is intrinsically non-axisymmetric and cannot be captured by 
our code. Additionally, several CTTSs are known where the higher 
order components are comparable to or stronger than the dipole 
components (e.g. BP Tauri, V2129 Oph, |Donati et al.|2008||20Ti] >. 
While the dipole may dominate the dynamics at larger distances in 
the disc, these higher order components play a major role in the dy- 
namics near the star. Additionally, axisymmetric simulations do not 
capture cases such as the one where the dipole field on the western 
hemisphere of the star is deflected in the opposite direction than 
the field on the eastern hemisphere. However, all of these effects 
would likely enhance the one-sided nature of the outflow as they 
contribute to the asymmetry of the magnetosphere. Further MHD 
simulations in full 3D are necessary to study these magnetospheric 
dynamics. 

Weak propellers: It is possible for an accreting star to be in the 
weak propeller regime of accretion where the star rotates faster than 
the inner disc, but not fast enough to drive an outflow. This occurs if 
the centrifugal force at the disc-magnetosphere boundary is insuf- 
ficient to overcome both the star's gravity and the magnetic tension 
opposing field line inflation. In this case some matter may still be 
ejected, but then it falls back onto the disc at larger distances from 
the star instead of flowing out of the system. In this sense, a weakly 
propelling star may continuously recycle the inner disc, causing the 
disc to evolve into a trapped state in which the accretion onto the 
star occurs in rare episodic bursts (D'Angelo & Spruit 2012). In 
the weak propeller regime, the disc may also form transient funnel 
flows and the outflow may be launched as magnetospheric ejections 
(Zanni & Ferreir af2013| l. However, if the propeller effect is espe- 
cially weak, no matter may be ejected at all (see also Romanova 
|et al.|2004| >. 

4.3 The effect of the MRI 

In prior works (e.g. Roman ova et al.|2 005 . Ustyugova et al.|2006| 
Romanova et al. 20091, propeller-driven outflows were observed 
in an a-disc where both the viscosity and diffusivity were mod- 
eled with numerical terms proportional to the Shakura & Sunyaev 
( 1973 ) a parameter. In these prior simulations, relatively large vis- 
cosity coefficients of a v = 0.1-0.3 were taken in the main models. 
Additionally, a large a diffusivity of ad = 0.1 was taken in order to 
facilitate the penetration of the disc matter into the magnetosphere. 

In the current work, we study the propeller outflow phe- 
nomenon in the case of a more physical MRI-driven disc which 
accretes due to the MRI-driven angular momentum transport. In 
contrast to the previous studies, no viscosity or diffusivity terms 
were added to the MHD equations, resulting in a more realistic 
model of the disc accretion. However, the MRI does not provide 
a source of diffusivity near the magnetically-dominated magneto- 
sphere and the inner disc matter only penetrates the magnetospheric 
field lines due to the very small numerical diffusivity in our MHD 
code. Nonetheless, our simulations show that the propeller-driven 
outflows persist even with the small effective diffusivity. In actual- 
ity, the diffusivity may be much larger in nature due to 3D instabil- 
ities such as the Rayleigh-Taylor and Kelvin-Helmholtz instability 
which act at the disc-magnetosphere boundaries (Romanova et al. 
[20081 [KuTkarni & Romanova|2008> . 

In Appendix [A] we study the effects of a larger diffusivity 
on the propeller outflow by applying an a-diffusivity to the matter 



near the disc-magnetosphere boundary. The larger diffusivity gen- 
erally increases the outflow rate, but only by a factor of two even in 
the most diffusive cases. The dynamics of the disc-magnetosphere 
boundary is poorly understood and hence, the effective diffusivity 
in the inner-disc is also poorly known; however, the simulations 
show that the diffusivity does not strongly affect the physics of the 
outflow. 



5 CONCLUSIONS 

In this work, we have used high-resolution simulations to study 
propeller-driven outflows from a realistic disc which accretes due to 
MRI driven turbulence. In the propeller regime of accretion, the star 
(and its magnetosphere) rotates faster than the Keplerian angular 
velocity of the inner accretion disc, resulting in a large centrifugal 
barrier at the disc-magnetosphere boundary. The main conclusions 
of this work are as follows: 

(i) Stars accreting in the propeller regime can launch centrifu- 
gally driven outflows. The rapidly-rotating magnetosphere creates 
a large centrifugal barrier which blocks matter from accreting onto 
the star. Instead, the accreting matter accumulates at the disc- 
magnetosphere boundary over time. Some of this matter diffuses 
through the outer magnetosphere into the more rapidly rotating in- 
ner regions where it picks up angular momentum from the field 
lines. If the centrifugal force is sufficiently large, the matter be- 
comes super-Keplerian and inflates the stellar field, creating a path 
for the matter to flow out of the system. As the matter flows away 
some of the field lines reconnect, releasing magnetic flux and mat- 
ter from the magnetosphere and truncating the outflow into discrete 
plasmoids. This matter dominated propeller wind is ejected with an 
half-opening angle of 45° and gradually collimated at larger dis- 
tances. 

The remainder of the matter which does not diffuse into the mag- 
netosphere instead undergoes an episodic accretion-ejection cycle. 
As the matter accumulates at the disc-magnetosphere boundary, the 
magnetosphere is slowly compressed towards the star. Eventually, 
the gravitational force becomes dominant and the matter deflects 
the magnetosphere to one side and accretes onto the star as a mag- 
netospheric funnel flow. A significant fraction of this accreting mat- 
ter is driven into an outflow, resulting in an enhanced outflow rate. 
Once the accumulated matter is depleted, the stellar magnetosphere 
quickly re-expands to halt the accretion onto the star. The enhanced 
outflow results in the ejection of a larger amount of magnetic flux 
which helps to collimate the matter, giving a narrower half-opening 
angle (compared to the non-accreting phase) of 20-40° at the edge 
of the simulation region. 

(ii) The propeller outflow has two main components: a high den- 
sity, medium velocity, matter dominated wind and a low density, 
high velocity, magnetically dominated Poynting jet. The Poynt- 
ing jet is magnetically accelerated and collimated by the helically 
wound magnetic field lines connecting the star to the corona. It 
carries substantial amounts of magnetic energy and angular mo- 
mentum away from the star and has a half-opening angle of 5-15° 
at the edge of the simulation region. The slower moving, matter 
dominated propeller wind is launched from the disc-magnetosphere 
boundary by the centrifugal force and helps to remove matter (and 
angular momentum) from the disc. Despite being matter domi- 
nated, the propeller wind also carries a substantial amount of mag- 
netic flux which acts to weakly collimate the outflow. We ob- 
serve weak collimation in our simulation region, but the magnetic 
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hoop-stress will continue to collimate the flow resulting in a well- 
collimated jet at large distances from the star. 

(iii) Stars in the propeller regime experience a spin-down torque 
through the magnetic field lines connecting the star to the slower 
rotating disc and corona and lose energy through the Poynting jet 
and propeller wind. The more rapidly rotating stars have higher 
spin-down rates and in all three of our models, we measure charac- 
teristic spin-down timescales of around 1 Myr for a typical CTTS 
with a strong 1 kG surface field (other parameters can be found in 
Table [Tj. 
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APPENDIX A: DIFFUSIVITY AND OUTFLOWS 

The MRI turbulence provides a mechanism for angular momentum 
transport in the disc and permits inward accretion from larger dis- 
tances. However, it does not act as a source of diffusivity at the 
disc-magnetosphere boundary (Romanova et al.|20lT| l. In spite of 
the MRI-driven turbulence, the accretion disc is matter-dominated 
and its dynamics are similar to a non-magnetized disc. The inner 
disc rotates at the Keplerian velocity whereas the magnetosphere is 
super-Keplerian and they do not interact unless there is a mecha- 
nism for diffusivity which provides penetration of the disc matter 
through the magnetosphere. 

In the ideal MHD simulations presented in the main text of 
this work, the matter in the inner disc gradually penetrates through 
the external magnetosphere due to numerical diffusivity. Fig. |A1| 
shows several snapshots of the accumulation of matter in the inner 
disc, the slow diffusive penetration through the external field lines 
of the magnetosphere, and the resulting outflow launched along the 
inflated stellar field lines. In our ideal simulations, the numerical 
diffusivity is low because the grid resolution is very high near the 
disk-magnetosphere boundary. In real systems, diffusive penetra- 
tion is determined by the magnetic Reynolds number: 

i?Ar 

Rm = , (Al) 

where v is the radial velocity of the inner disc matter, Ar is 
the characteristic size of compressed layer of the magnetosphere 
through which the inner disc matter should diffuse, and r/ m is the 
coefficient of the magnetic diffusivity. At the disc-magnetosphere 
boundary of the propelling star, the inner disc is almost stopped 
by the magnetosphere and therefore the radial inward velocity v 
is small. Additionally, the magnetosphere is compressed by the in- 
ner disc meaning that the characteristic scale Ar is also small, in 
particular during episodes of matter accumulation. Therefore, the 
condition Re m < 1, which is necessary for diffusive penetration, 
is satisfied even at small r\ m as we observe in our simulations. 

The diffusivity of the stellar magnetosphere is one of the most 
poorly understood properties of young stars and may be larger than 
the numerical diffusivity present in the code. One of the main lim- 
itations of our axisymmetric simulations is that we cannot directly 
simulate mechanisms such as the Rayleigh-Taylor and Kelvin- 
Helmholtz instabilities which may be responsible for the mixing 
of the matter and field at the disc-magnetosphere boundary l |Arons| 
|& Lea|1976| >. In 3D, these mechanisms may serve as the primary 
source of diffusivity in the magnetosphere. The Kelvin-Helmholtz 
instability may be important because the inner disc matter rotates 
slower that the external magnetosphere. In the propeller regime, the 
Rayleigh-Taylor instability cannot act in its usual form because the 
presence of an inward pointing effective gravity is required ( Spruit 
et al. 1995. Kulkarni & Romanova 2008). However, it may operate 
in the opposite direction, driven by the outward pointing centrifu- 
gal force. The enhancement of diffusivity can strongly influence the 
rate of accretion and outflows in the propeller regime. 

To investigate the role of diffusivity in our axisymmetric sim- 
ulations, we developed a non-ideal version of our MHD code 
by adding diffusivity terms to our MHD equations. The diffusiv- 
ity coefficient is proportional to the Shakura & Sunyaev| ^1973} 
a— parameter. We take a — ay inside the radius r < 5 where the 
inner disc interacts with the magnetosphere, and a = at larger 
radii so that the MRI-driven accretion proceeds normally in the rest 
of the disc (high diffusivity damps the MRI). We also introduce 
the threshold density pd below which the diffusivity is not applied. 
Therefore, we take into account only the higher-density matter of 



the inner disc and exclude diffusive penetration of the coronal mat- 
ter. This restriction helps to imitate the Rayleigh-Taylor instability 
at the disc-magnetosphere boundary which develops more easily at 
sufficiently high densities (this is because low-density small-scale 
filaments coalesce forming larger and denser filaments which pene- 
trate through the magnetosphere, e.g. Kulkarni & Romanova 2008) 

As a base, we take the case of the propelling star with coro- 
tation radius r cor — 1.5 and vary the diffusivity ay as well as 
the matter density above which the diffusivity operates pd- In the 
first series of runs, we choose a relatively small threshold density 
pd — 0.01 and perform simulations at three different diffusivities: 
ad — 0.01, 0.1 and 1 (see set 1 in Table | A l) . In another series of 
runs, we take large diffusivities: ad — 0.1 or 1 and high threshold 
densities: pd = 1 and pd = 5 (set 2 in Table |Al| >. This second set 
of simulations is meant to study the Rayleigh-Taylor mechanism in 
the magnetosphere. 

For each simulation, we calculate fluxes of matter onto the star 
M* and into the wind M w i n< j. The fluxes into the wind are calcu- 
lated through the surface z = 20, r — 20 for velocities v p > 0.2. 
The top panels of Fig. |A1| show that at small threshold densities 
Pd — 0.01, the flux onto the star systematically decreases with 
increasing oj and becomes less spiky, whereas the matter flux to 
the wind becomes somewhat larger. This dependence of fluxes on 
diffusivity can be easily understood: at larger values of ay, mat- 
ter diffuses more efficiently through the field lines and is ejected 
faster than it can accumulate at the disc-magnetosphere boundary. 
In cases of low diffusivity, the situation is opposite: the accreting 
matter does not efficiently diffuse through the field lines, instead 
accumulating at the disc-magnetosphere boundary where it eventu- 
ally overflows onto the star. This episodic accumulation-overflow 
results in the spikes in accretion which are observed. 

We calculate averaged fluxes onto the star and into the wind 
using Equation [T] We expect that a higher diffusivities will re- 
sult in a higher outflow efficiency and indeed, we find that the 
accretion rate onto the star (M*) systematically decreases with 
diffusivity in the first set of simulations. Matter flux to the wind 
(M w i n d} increases for ad = 0.1 then decreases again when the 
diffusivity becomes too high at ad = 1. The angular momentum 
flux from the star, (L*) (which determines the spin-down of the 
star), correspondingly increases, then decreases. It has a maximum 
value of 2.7, which is only 1.7 times larger than that in case of 
ad — 0. Therefore, if additional diffusivity were to be included in 
our base models, the spin-down rate of the star may increase, but 
only by a factor of ~2. As expected, the efficiency of the propeller 
(Afwind)/((A^wind) + (M*)) systematically increases with ad- 

In our second set of simulations, we only enable the diffusivity 
if the local density is larger than a density threshold value, pd- In 
these cases, matter accumulates at the disc-magnetosphere bound- 
ary until it exceeds the threshold density of pd', then it efficiently 
diffuses into the magnetosphere, with some fraction accreting onto 
the star and the remainder ejected into the wind. The bottom pan- 
els of Fig. |A2| show that the accretion and a significant number 
of the ejections proceed in spikes, reflecting the long episodes 
of matter accumulation at the disc-magnetosphere boundary. The 
timescale between major spikes is determined by the timescale of 
the matter accumulation — typically on the order of tens of dynam- 
ical timescales. This timescale depends on both of the parameters, 
ad andpd- 

Analysis of this section shows that the larger diffusivities ex- 
pected in more realistic three-dimensional simulations may lead to 
somewhat stronger outflows and shorter spin-down timescales of 
the propelling stars. These larger diffusivities model the mixing in- 
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Figure Al. Diffusive penetration through the magnetosphere. The snapshots show the gradual diffusive penetration of the disc matter through several field 
lines of the outer magnetosphere. The color background shows the logarithmic density distribution with red indicating the maximum density and light-cyan 
indicating the smallest density in the disc. There is also an even lower density corona which is shown in white. 





base case 


set 1 : low p,; 


set 2: high p^ 


ad 

Pd 

(M wind ) 

(M,) 

<A'/ Wi „d> 


0.0 

0.0 

0.27 

0.09 

0.74 

1.7 


0.01 0.1 1.0 
0.01 0.01 0.01 
0.26 0.4 0.34 
0.11 0.08 0.05 
0.70 0.83 0.86 
2.0 2.7 2.25 


0.1 0.1 1 
1 5 5 
0.34 0.16 0.24 
0.12 0.08 0.07 
0.74 0.67 0.77 
2.0 1.70 1.75 


<M wind > + <M„> 

CM 



Table Al. Summary of results for the diffusive runs. All runs were performed for r COI = 1.5. ay is the diffusivity parameter, p^ is the threshold density for 
diffusivity, (Myriad) an d (M*) are time-averaged values of matter flux to the wind and to the star, respectively. The ratio (JW w i Iu j)/((-^wiiid) + (M*)) is 
the outflow ejection efficiency of the propeller, (i* } is the angular momentum flux from the surface of the star. 



stabilities at the disc-magnetosphere boundary which we do not see 
in our 2.5D simulations. Our investigation shows that the difference 
between the cases with very low diffusivity and high diffusivity is 
not substantial, differing only by a factor of 2. The reason is that 
in both cases (for the same accretion rate in the disc) roughly the 
same amount of matter is ejected to winds due to propeller mecha- 
nism. In the less diffusive runs, the ejection occurs episodically in 
bursts whereas in the more diffusive runs the process is smoother. 
On average, however, the matter fluxes are comparable. 



APPENDIX B: STRESSES IN THE MRI DISC 

The matter accretion onto a star in the propeller regime is deter- 
mined by the turbulent stress in the disc. The integrated matter 
(subscript "m") and magnetic field (subscript "/") stresses in the 
co-moving frame are given by the vertically averaged angular mo- 
mentum fluxes (Equation [5J: 

( T ™) = T^E I dzpv r v^ - (pv r )(v^) (Bl) 



and 

(Tf) = 
where 



2H j 

1 
~2H 



dz 



B r B$ 

4.7T 



(Vfi) = — / dzpv^, (pV r ) 



1 

2H 



dzpv r , 



(B2) 



(B3) 



E = J dzp is the disc surface density and 2H is the total thickness 
of the disc which we set using a density threshold. The matter and 
magnetic pressures are given by 



(Pm) = 



2H 



dzP, (P f ) = 



1 
2H 



dz^^. 
8tt 



(B4) 
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Figure Bl. Radial distribution of stresses in the disc. The left panel 
shows the z-averaged stresses and pressures inside the disc. Outside the 
magnetospheric region, the magnetic stress dominates over the matter 
stress. The right panel shows the computed a-parameters. 



The standard a— parameters associated with matter and magnetic 
stresses are defined as the height-averaged stress divided by the 
height-averaged pressure 



(Tm) 
(Pn,)' 



a f 



(Tf) 



Fig. |B1| shows the radial distribution of stresses, pressure and 
a— parameters for the simulation. The magnetic stress Tf roughly 
an order of magnitude larger than the matter stress throughout the 
disc. Near the magnetosphere, the magnetic pressure dominates 
over the matter pressure, but is smaller in the external parts of the 
disc. Correspondingly, the magnetic a-parameter aj is larger than 
the standard a-parameter determined by Reynolds stress, a m . In 
the magnetospheric region (r < 8), both a-parameters are very 
large because both stresses are large. Further away from the star, 
we find a f fa 0.04 - 0.1 and a m » 0.001 - 0.08 indicating 
that the radial angular momentum transport and inward accretion is 
determined by the magnetic stress. 
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Figure A2. Matter fluxes to the star (top curves) and to the wind (bottom curves), calculated through the surface r = 20, z = 20 and restricted by velocity 
v > 0.2, for different values of ay and density thresholds p^ (at which the diffusivity is applied). 
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